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ABSTRACT 

UV radiation from massive stars is thought to be the dominant heating mech- 
anism of the nuclear ISM in the late stages of evolution of starburst galaxies, cre- 
ating large photodissociation regions (PDRs) and driving a very specific chem- 
istry. We report the first detection of PDR molecular tracers, namely HOC" 1 ", 
and CO + , and confirm the detection of the also PDR tracer HCO towards the 
starburst galaxy NGC 253, claimed to be mainly dominated by shock heating and 
in an earlier stage of evolution than M82, the prototypical extragalactic PDR. 
Our CO + detection suffers from significant blending to a group of transitions of 
13 CH 3 OH, tentatively detected for the first time in the extragalactic interstellar 
medium. These species are efficiently formed in the highly UV irradiated outer 
layers of molecular clouds, as observed in the late stage nuclear starburst in M 82. 
The molecular abundance ratios we derive for these molecules are very similar 
to those found in M82. This strongly supports the idea that these molecules 
are tracing the PDR component associated with the starburst in the nuclear re- 
gion of NGC 253. The presence of large abundances of PDR molecules in the 
ISM of NGC 253, which is dominated by shock chemistry, clearly illustrates the 
potential of chemical complexity studies to establish the evolutionary state of 
starbursts in galaxies. A comparison with the predictions of chemical models for 
PDRs shows that the observed molecular ratios are tracing the outer layers of 
UV illuminated clouds up to two magnitudes of visual extinction. We combine 
the column densities of PDR tracers reported in this paper with those of easily 
photodissociated species, such as HNCO, to derive the fraction of material in 
the well shielded core relative to the UV pervaded envelopes. Chemical models, 
which include grain formation and photodissociation of HNCO, support the sce- 
nario of a photo-dominated chemistry as an explanation to the abundances of the 
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observed species. From this comparison we conclude that the molecular clouds in 
NGC253 are more massive and with larger column densities than those in M82, 
as expected from the evolutionary stage of the starbursts in both galaxies. 

Subject headings: galaxies: abundances — galaxies: ISM — galaxies: starburst - 
galaxies: individual(NGC 253) 
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Introduction 



Intense UV radiation from massive stars is one of the main mechanisms responsible 
for the heating of the interestelar medium in the nuclear region of starburst galaxies. This 
mechanism is particularly important in the latest stages of starburst (SB) galaxies where 
the newly formed massive star clusters are responsible for creating large photodissociation 
regions (PDRs). This is the case for the prototypical SB galaxy M82, where the large 
observed abundances of molecular species such as HCO, HOC + , CO + , and H 3 + are 
claimed to be probes of the high ionizatio n rates in large PDRs formed as a consequenc e 



of its extended evolved nuclear starburst (iGarcia-Burillo et al 



van der Tak et al. 



2008|). 



2002 



Fuente et al 



2006 



Observational evidences point to a significant enhancement in the abundance of HOC 4 " 
in regions with large ionization fractions. The abundance r atio [HCO + |/[HOC" f |= 270 is 



found in the prototypical Galactic PDRs of the Orion Bar (I Apponi et al. 



or ey en lower abundance ratios are observe d in the PDRs NGC7023 (50-120, 



2003|), SgrB2(OH) and NG C2024 (360-900, 



and the Horsehe ad (75-200 



Ziurvs fc Apponi 



Goicoechea et al 



1995 



1999h. Similar 



Fuente et a 



Apponi fc Ziurvs 



19970 



Liszt et al. 



2009|), as well as in diffuse clouds (70-120, 



20041 ). This is in contrast with the much larger ratios of 3> 1000 found in dense 



molecular clouds well shielded from the UV radiation. However, these low HCO + /HOC 



;arge values o 



this ratio of > 2000 are found 



1999 



Savage fc Ziurvs 



200J). 



ratios are not found in other galactic PDRs 
in the PDRs M17-SW, S140, and NGC2023 flApponi et al. 
The HCO molecule has also been observed to be a particularly good tracer of the PDR 
interfa ces. Low ratios of [HCO + ]/[HCO]~ 2.5 — 3 are found in prototypical Galactic 



PDRs (ISchenewerk et al. 



1988 



Schilke et al. 



200 lh . The large HCO abundance (> 10 



altogether with the low ratio [HCO + ]/[HCO]~ 1 in the Horse head PDR is claim ed to be 



a diagnostic for an ongoing FUV-dominated photochemistry (jGerin et al. 



2009|). CO+ is 
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also claimed to be particularly prominent in the chemical modeling of PDRs and high 



abundances of this mo 



([Sternberg fc Dalgarno 



ecule appear to be correlate d to similar enhancements of HOC 



1995 ; 



Savage fc Ziurvs 



1-10 are observed in a number of PDRs ( 



the Horsehead PDR (jGoicoechea et al 



2004). rCO+l/fH OC+l ratios in the range of 



Savage fc Ziurvs 



2009|) 



200J), but only of > 0.1. towards 



As mentioned above, this set of PDR probes has been extensively studied towards 
M 82. However, no such complete studies have been carried out toward s other prototypic al 



galaxies, but for the de tection of HCO and HO C+ towards NGC 1068 (lUsero et al.ll2004l ) 



and H 3 0+ in Arp220 (; 



van der Tak et al 



2008T ). M82 and NGC 253 are the brightest 



lumin osities and star formation rates of about ~ 3 M yr 1 (jOtt et al 



showing very similar IR 


(Ott et al. 


2005; 


Minh et al. 



20071 ). However, both galaxies show very different chemical composition. The chemistry and 



to a large extend the heating in th e central region of N GC 253 is believed to be dominated 



by large scale low velocity shocks fjMartm et al 



2006b! ) . The similar chemical composition 



found in the nuclear region of NGC 253 to that in Galactic star forming molecular complexes 



Martin et al. 


2003. 


2005. 


2006bl 



Furthermore, our recent observations of the PDR component as traced by the e asily 



2008|) 



photodissociated HNCO molecule towards a sample of galaxies (iMartin et al. 
showed the non-detection of HNCO in M82, at a very low abundance limit. This low 
HNCO abundance supports the scenario that the PDR chemistry dominates the molecular 
composition of the ISM in this galaxy. However, from the HNCO measured abundance in 
NGC 253, it would be placed in an intermediate stage of evolution where photodissociation 
should be starting to play a significant role in driving a UV-dominated chemistry which has 
not been yet identified towards this galaxy. The presence of a significant PDR component 
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in NGC 253 claimed from the HNCO abundance is also inferred from the simi 



of th e atomic fine stru cture line intensities from PDR tracers like CII and 01 (ICarral et al- 



ar intensity 



1994 



Lord et al 



1996|) observed in both M82 and NGC 253. 



In this paper we present the first detection of PDR molecular tracers HO C 4 " and CO 



Sage fc Ziurys 



19951 ) in the central 



and confirm the detection HCO (tentatively detected by 
region of NGC 253 which allows the evaluation of the influence of the photodissociation 
radiation in the nuclear ISM of this SB galaxy. The results presented here support the 
scenario of the presence of a significant PDR component and clearly show the potential of 
molecular complexity in estimating the contribution of the different heating mechanisms of 
the ISM in the nuclei of galaxies. 



2. Observations and Results 

The observations presented in this paper were carried out at the IRAM 30 m and 
JCMT telescopes on Pico Veleta, Spain, and Mauna Kea, USA, respectively. 



2.1. IRAM 30 m 

The IRAM 30 m observations were performed in symmetrical wobbler switched mode 
with a frequency of 0.5 Hz and a beam throw of 4' in azimuth. The 516 x 1 MHz filter banks 
were used as spectrometers. 

We have observed the transitions of C 18 J = 1 - (109.782 GHz), HCO+ J = 1 - 
(89.188GHz), HOC+ J = 1 - (89.487GHz) and 3 - 2 (268.451GHz), and the HCO 
lo,i — Oo,o (86.670 GHz). Beam sizes at these frequencies were 22", 28", and 9". The nominal 
position for the observation was aj2ooo = 00 h 47 m 33 f !3, 5j2ooo = — 25°17'23" for HCO 4 " and 
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HOC + , matching up the position used for the 2 mm line survey (IMartm et al.ll2006bl ). The 
data on C 18 and HCO were centered at a J20 oo = 00 h 47 m 33?5, 5 J20 oo = -25°17'27". The 
two positions are separated by < 5" which, considering the 3 mm beam sizes, should have a 
negligible influence in the relative intensities. Double Gaussian profiles have been fitted to 
all observed transitions and the corresponding derived fitting parameters are summarized 
in Table [TJ 

Fig. Q] shows the simultaneously observed J = 1 — features of HCO + and HOC + 
compared to the C 18 J = 1 — line profile. Although at this frequency the SIS receivers 
image band rejection is larger than 20 db, we observed the HOC" 1 " line tuned to two 
different velocities (250 and 500 kms -1 ) in order to confirm that the observed profile was 
not line emission coming from the upper side band. Fig. [T] shows the average of both 
observations. The HCO + J = 1 — was detected at the edge of the band covered in the 
500 kms -1 tuning of HOC + . Although pointing accuracy was of the order of 3", the 
different shapes observed between the HCO + and HOC + is attributed to a small change in 
the pointing position during the two observations. This effect accounts for an uncertainty 
in the integrated intensity of < 10%. Even though the HCO + feature was not completely 
covered by the backen d, the integrated inte nsity we derive is in agreement within 5% with 



that of theobserved by 



Nguyen et al. 



(119921 ) and we observe the line shape to be consistent 
with that of C 18 0. As indicated in Tabled! only a coarse upper limit to the detection of 
HOC + J = 2 — 1 was obtained. 

The HCO J = 1 — emission was observed in the same window as SiO 2 — 1 and 
H 13 CO + 1—0 and appears slightly blended to the latter. With a significantly improved 
signal-to-noi se ratio, we co n firm t he previous tentative detection of this HCO transition 



reported by 



Sage & Ziurvsl (119951 ) with the NRAO 12 m 



main beam brightness temperature from 



teles cope. Moreover, using the 



Sage Sz Ziurvsl (119951 ) of ~ 1 mK with at 72" beam 
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and our observed ~ 4 mK with a 28" beam we can make an estimate of the emitting source 
extent of > 20". The double Gaussian profiles fitted to each species were constrained 
to have similar linewidths. The resulting fitted line positions agree within the errors to 
those expected from the rest frequencies of each line. Fig. [2] shows the results of the fit 
superimposed on the observations as well as the position of the hyperfine structure lines 
of HCO. Only the brightest of the group (F — 2 — 1) has been taken into account for the 
fit. Assuming optically thin emission, the F = 1 — and F — 1 — 1 transitions (at 86.708 
and 86.777 GHz) are expected to show an intensity half of the main transition but they are 
completely blended to the H 13 CO + emission. The F = — 1 transition at 86.805 GHz is 
expected to be even fainter by a factor of 5, well below our detection limit. Fig. [2] shows 
in dotted line a synthetic spectrum of HCO assuming one velocity component centered at 
255km s" 1 with a linewidth of 192km s _1 (as derived if only one component is fitted to the 
spectrum from the other lines) and a peak intensity of the HCO F = 2 — 1 line of 3.6 mK. 
This shows that the fainter HCO hyperfine transitions may account for up to a 10 — 20% of 
the H 13 CO + integrated intensity. 



2.2. JCMT 



JCMT observations were performed in beam switched mode with a frequency of 1 Hz 
and beam throw of 2' in azimuth. The ACSIS digital autocorrelator spectrometer was used 
with a bandwidth of 1600 MHz providing a resolution of ~ 1 MHz. 

We have used the receiver A3 to observe the CO + transition at 236.062 GHz. At 
this frequency, the beam size of the telescope is 21" and the main beam efficiency 0.69. 
The observations were carried towards the nominal position «j2ooo = 00 h 47 m 33 f !l, S d 



-25°17'18" (radio continuum position, 



Douglas et al. 



; J2000 



19961 ). As seen in the HC 3 N 



J = 25 — 24 profile, most of the emission is observed from one of the velocity components 
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at this position, which is due to the JCMT observed position being ~ 6" and 10" away 
from those observed with the IRAM 30 m, respectively. This position is half beam away 
from the positions observed with the IRAM 30m, so the abundance ratios derived from this 
observation might be affected by a larger uncertainty of up to a factor of 2. However, this 
effect might be attenuated by the emission being extended over scales of > 20". 

The CO+ 5/2 - 3/2 F = 2 - 1 and 3/2 - 1/2 F = 2 - 1 transitions are clearly detected 
above the noise level (~ 1.5 mk in 30kms _1 channels). However, we observe its profile 
significantly blended to that of the group of transitions of 13 CH 3 OH J = 5 — 4. This 
overlap was not a problem in the case of the CO + detection towards M 82 due to the 
significantly lower abundance of CH3OH towards this galaxy (IMartm et al.ll2006al ). The 
CO + 3/2 — 3/2 F = 2 — 1 component was not detected due to its low relative intensity. 
We have fitted the observed profile with a single Gaussian component CO + synthetic 
spectra with radial velocity and linewidth fixed from those derived from HC3N. The relative 
intensities of the CO + components were also fixed to those expected from optically thin 
emission under under local thermodynamic equilibrium conditions. The derived line profiles 
parameters are presented in Table [TJ Additionally, we simultaneously fitted a synthetic 
13 CH30H spectra, using the 18 13 CH30H transitions in the observed frequency range, using 
the same constraints as for the CO + line fit. The fit reproduces most of the observed 
features and shows that the emission from 13 CH 3 OH may explain the observed non Gaussian 
CO + profiles. Only the parameters derived for the three most intense components in the 
group are given in Table [TJ This is the first time that the 13 C isotopologue of methanol is 
detected towards an extragalactic source. Regarding the accuracy of the fitted parameters 
presented in Table [TJ the integrated line intensities derived for CO + and 13 CH 3 OH are 
likely underestimated by ~ 20% due to the baseline determination. In the next Section, we 
discuss the detection of 13 CH30H in the context of the derived abundances with respect to 
those of the main methanol isotopologue. 
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Molecular abundances and ratios 



We have estimated the fractional abundances of the newly observed species in NGC 253 
assuming optically thin emission, LTE conditions, and similar spatial distribution for all 
species. Under these assumptions, we have calculated the column densities of H 13 CO + , 
HOC + , HCO, and CO + for an excitation temperature T ex = 15 ± 5K and an estimated 
source extent for each velocity component of 10". The T ex = 15 ± 5K is assumed based 
on the ave rage rotational temp eratures derived from most of the species detected towards 



NGC 253 flMartm et al 



2006bl ). Indeed the non detection of HOC + 3 — 2 implies low 
excitation temperatures of T cx ~ 10 K. Both the excitation temperature and the emission 
extent have an important impact in the absolute derived column densities by up to a factor 
of 2, however, the fractional abundances and abundance ratios are mostly independent of 
these assumptions. We assume that the emission extent is similar for all observed species. 
Table [2] presents the column densities and fractional abundance ratios with respect to H 2 for 
all the species. The total H2 column density has been derived from the C 18 Q column densit y 



for ea ch velocity component assuming an isotopic ratio of 16 0/ 18 = 150 ([Harrison et al 
19991 ) and a CO/H2=10 -4 . The relative abundances derived for all molecules are presented 



in Table [2j Additionally, abundance ratio of H 13 CO + with respect all species is also 
presented in Table [2j The errors in the derived column densities take into account the 
statistical error of integrated intensities and the uncertainty in the assumed excitation 
temperature. These errors are subsequently propagated to the abundance ratios. The 
HCO + abundance has been derived from that of H 13 CO + to avoid the opacity effects 
affecting the m ain isotopologue. Indeed, if the 12 C/ 13 C ratio of ~ 40 derived for NGC 253 



and NGC 4945 (IHenkel et al. 



1993 



1994 



Martin et al. 



20051 ) applies to all these galaxies, 
an average opacity of THnn+ J=i-n > 1 is derived from the HCO + /H 13 CO + J = 1 — line 



ratio (this work and 



Usero et al. 



20041 ). As a consequence, ratios may be underestimated by 



a factor of ~ 2 if derived from the main isotopologue. The CO + column density presents 
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a significant uncertainty due to its blending with the newly detected isotopologue of 
methanol, 13 CH 3 OH. However, the CO + abundance should not be affected by more than a 
factor of 2. 

We have derived a column density for 13 CH30H of ~ 1.6(1.0) x 10 u cm" 2 . If w e 



compa re this column density with that of the main isotopologue from (IMartm et al. 
2006bl ) it res ults in a CHaOH/^ CHaOH ratio of 12 ± 7, significantly lower than the ratio 



12 C/ 13 C=40 flHenkel et al 



19931 ). Both the difference of 6" in the observed positions and 



a different filling factor of CH3OH might account for part of this difference. However, 
th e integrated intensitie s me asured for the meth a nol gr oup of transitions at 145.1GHz 



by 



Martin et al. 



fl2006bh and 



Hiittemeister et al. 



(119971 ) at positions differing by > 13" 



show a variation of < 6% so the difference in positions is not likely to contribute to 
this differen ce. Thus, opacity is likely t he dominant effect as observed in the Galactic 



center (GC, 



Requena- Torres et al 



2006|). From the 12 C/ 13 C ratio in NGC253 we derive 



a fractional abundance of methanol o 



Galactic center (IRequena- Torres et al. 



~10 7 , close to the abundances observed in the 
20061 ). Methanol is, after C O and NO, the most 



abundant molecule in the nucleus of NGC253 (IMartm et al. 



2006H ). 



4. Discussion: The PDR component in NGC 253 

4.1. NGC 253 in context 

Table m shows the HCO+/HOC+, HCO+/HCO, and HCO+/CO+ abundance ratios 
resulting from our measurements in NGC 253 compared to those of the similar SB galaxy 
M82, and the Seyfert 2 with nuclear SBs, NGC 1068, and NGC 4945, together with 
prototypical galactic PDRs, where observations of these species have been reported. For 
the shake of consistency, the abundance ratios of the other galaxies have been calculated 
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from the available line profiles obtained from the literature. As already explained before, 
for this comparison we used H 13 CO + to derive the HCO + column densities i n order to 



avoid the opacity effects. We have assumed the isotopic ratio of C/ C=40 (IHenkel et al 



19931 ) to derive the column densities of the main isotopologue. For M82, the t 



and H CQ + /CO + abundance ratio have been derived from the observations by 



CO+/HOC+ 



Fuente et al. 



( 120061 ) towards the Eastern molecular lobe. Our measured ratios are ~ 59 and ~ 0.8, which 



are > 30% larger that the ratios derived by 



Fuente et al 



(120061 ). T his difference is due 



to th e significant missing flux of the H 13 CO + interferometric maps (jGarcia-Burillo et al 
20021 ) they used for comparison, as well the higher 12 C/ 13 C ratio they used to compare 



with single dish HCO + observations. By c omparing the convolved in tegrated intensity 



from the H 13 CO + i nterfe rometric maps by (jGarcia-Burillo et al 



lobe 



Fuente et al. 



(120061 ) and the single dish data of by 



2002) towards the western 



Mauersberger fc Henkell (ll99ll ) 



towards a nearby posit ion, we have estimated a > 50 % missing flux. Thus we used the 



H 13 CO+ 1 - data by 



Mauersberger fc Henkell (ll99ll ) in our measured ratios. The ratios 



towards NGC 106 8 were derived for the regions within the g alaxy where li nes in tensities 



were tabulated by 



Usero et al. 



(120041 ) . The observations from 



Wang et al 



(120041 ) were used 



to derive the HCO+/HCO ratio towards NGC 4945. 



4-1.1. PDR abundance ratios in SB galaxies 

We find that the three derived HCO+/HOC+, HCO+/HCO, and HCO+/CO+ 
abundance ratios in the two starbursts, NGC 253 and M82, are equivalent within the 
measurement errors. The ratios are also in reasonable good agreement to those found in 
galactic sources with similar FUV fluxes (see Table [3]) • Such high abundances ratios of 
HOC + , HCO, and CO + relative to HCO + have been claimed to be the evidence of M82 
being mostly dominated by photodissociation. We notice that the average HCO + /HCO 



13 



ratio in NGC253 are even lower than that measured in M82. In the case of HCO, the 
interferometric maps of M 82 clearly resolve the spatial variations in this ratio across the 
galaxy nuclear region. However, towards the region of peak HCO emission in the M 82 
maps, we find a ratio of H 13 CO + /HCO ~ 0.12 ± 0.04, equivalent to the average observed 
towards NGC253. Our data show that the ISM in the nuclear region in NGC253 must be 
significantly pervaded by a strong UV radiation flux from the massive star clusters formed 
i n the starburst, as also suggested by the study of the abundances of the HNCO / CS ratio 



( IMartm et al. 



20091 ). Moreover, these new observations would imply that photodissociation 



plays a similar role in the ISM heating of both NGC253 and M82. 



4-1.2. PDR abundance ratios in AGN galaxies 

Both HCO+/HCO and HCO+/HOC+ abundance ratios in NGC 1068 are different by 
a factor 2-3 from those of NGC 253 and M82. Furthermore, HCO+/HCO is also found 
to be up to a factor of ~ 2 lower in the ring of star formation than towards the nuclear 
region. Like in NGC 253, these ratios are consistent with the decrease i n the abundance o f 



molecules such as HNCO from the nuclear region to the starburst ring (IMartm et al 



The tentative detection of HCO in the circunnuclear disk (CND) of NGC 1068 (jUsero et al. 



2009). 



20041 ). suggest a rough H 13 CO + to HCO line intensity ratio of ~ 2 — 3, which turns into an 
abundance ratio in the range of ~ 0.09 — 0.13, closer to the values derived in SB dominated 
galaxies. This value at the CND can be significantly biased by the emission from the star 
forming ring covered at half-power by the 28" beam. Therefore, it is clear that this ratio 
does not significantly decrease towards the nuclear AGN in this galaxy with the typical 
angular resolution of 20" — 30". On the other hand, the ratio HCO + /HOC + is only a factor 
of 2 lower in the starbursts galaxies than towards the nuclear AGN in NGC 1068. From 
these observations, it is unclear whether photodissociation does play a major role in the 
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AGN dominated center of NGC 1068. Unfortunately, no observations of HOC + are available 
towards the SF ring in this galaxy. 



Similar to NGC 1068, the obscured Seyfert 2 nucleus in NGC 494 5 is surrounded by a 



starburst ring more prominent than in NGC 1068 (IGenzel et al. 



1998|). The HCO+/HCO 



ratio found in NGC 4945 is even lower than that fou nd in the oth e r gala xies. However, the 



fit to these lines was claimed to be very uncertain by 



Wang et al 



f l2004h . 



4.2. Comparison to PDR chemical models 



We have compared our observed frac tional abundanc es and abundance ratios with 



those predicted by the UCL_PDR model (IBell et al 



20061 ). The UCL_PDR code is a time 



and depth dependent one dimensional PDR model that si multaneously so lve the chemistry, 



thermal balance and radiative transfer within a cloud (see 

\5 „™-3 



Bell et al. 



2006 



for more details). 

We have adopted a hydrogen density of 10 5 cm -3 , a radiation field Go ~ 5000 in units 
of Habing field, and a cosmic radiation rate of 10~ 16 . The high density of 1 5 cm -3 is 



2008 



2009 



; Aladro et 



derived from the multiline analysis of CS and HC3N (IBayet et al. 
al. In Prep). Our estimates of the cloud structure (see Sect. 14. 31) depends on the averaged 
radiation field which might be different in both galaxies. The averaged radiation fields in 
both NGC 253 and M82 have been inferred from the fine str ucture lines and they are of 



2 x 1 4 and 10 3 , respectively, with large errors of a factor of 2 fjCarral et al 



1994 



Lord et al. 



19961 ). This would imply that the PDR envelope should be larger in the clouds of NGC 253 
than in M82. Since we do not aim to quantitatively model the particular abundances 
measured in NGC 253 but to investigate the physical conditions that would give rise to the 
wealth of observed molecules in starburst, the value of Go = 5 x 10 3 used is a geometric 
mean value derived from the fine structure lines in both galaxies. Given that our estimates 
of the A v for the PDR are based on this geometric mean, the expected changes in A v for 
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the two galaxies would be just a factor of 1.6. 

We have ran two different models: Model A is a standard time dependent gas-phase 
PDR model where the initial composition is atomic; while Model B is computed using a 
coupled dense core-PDR model where the diffuse material, initially also purely atomic and 
gaseous, collapses to reach a final density of 10 5 cm -3 . During the collapse the gas depletes 
on the grains forming icy mantles which remain on the dust until irradiation from a UV 
field is switched on, evaporation occurs and the typical PDR chemistry takes place. In 
both models the temperature is calculated self-consistently at each depth and time step 
by thermal balance. Fig. H] shows the predicted abundances and abundance ratios as a 
function of visual extinction (A v ) for Model A (left panels) and B (right panels). HNCO 
and CH3OH results are only shown for Model B in Fig. HI 

Observed abundances of HCO + and HOC + towards NGC 253 (shown as horizontal 
lines in the key of Fig. H|) are well reproduced by the models for very low extinction 
of A v ~ 1—2. The HCO abundance observed is a factor 2 — 6 above the maximum 
predicted by the model B. It is important to take into account that while we assumed a 
ratio CO/H 2 =10~ 4 to calculate the fractional abundances, this ratio is not constant in the 
models and hardly ever reaches this value. On the other hand, the abundance ratios of 
the observed molecules, unaffected by the hydrogen determination uncertainty, agree well 
with the model predictions. The model shows that the abundances of CO + , follows the 



same pattern as HOC + . The c orrelation between these two molecules was a 



by previous theoretical studies ([Sternberg fc Dalgarno 



1995 



Savage fc Ziurys 



so pre dicted 



2004h . Thus, 



CO + measurement allows us to confirm the effect of photodissociation suggested by the 
large abundance of HOC + . 

The ions HOC + , CO + and HOC + are mostly formed at the edge of the cloud and 
while the two models predict similar abundances for these ions, it is worth noting that 
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the observed abundances of other species such as HNCO, CH3OH, or HOCO + can only be 
explained with Model B. However, HNCO and CH 3 OH only reach observable abundances 
for A v ~ 5 magnitudes. This implies that, while photodissociation does play an important 
role in the chemistry of NGC253, the molecular clouds affected by the UV radiation must 
contain a dense core well shielded from the UV radiation and rich with gas phase icy mantle 
molecules like HNCO and CH 3 OH. 



4.3. The molecular clouds in SB galaxies 



Martin et al.l (120081 ) have used a comparison with the GC molecular clouds to propose 
another PDR diagnostic based on the relative abundance of HNCO to CS. The large 
variation of the HNCO/CS abundance ratio between UV radiated clouds to those well 
shielded clouds only affected by shocks was interpreted as the fast photodissociation of the 
fragile molecule HNCO, efficiently produced on the icy mantles a nd delivered into ga s phase 



Martin et al. 



fl2009h found 



by low velocity shocks. Similarly, in a sample of nearby galaxies 
changes of nearly two orders of magnitude from the shock dominated chemistry in M 83 and 
IC342 to UV dominated chemistry in M82. The extremely low HNCO abundance in M82 
and the large abundances of HCO, HOC + and CO + support the idea that the HNCO/CS 
ratio is a measure of the relative importance of the UV heating to shock heating and the 
evolutionary state of the starburst in galaxies. 

The detection of HCO, HOC + and CO + in NGC253 with similar column densities and 
abundance to those in M82 suggest that the PDR component is similar in both galaxies 
as suggested by the similar atomic fine structure line intensities in both galaxies. The 
results of model B confirm the observational trends observed in HNCO and the other 
PDR molecules in galaxies. For molecular clouds with A v ~ 1 — 2 (i.e column densities 
of 1 — 2 x 10 21 cm -2 ) illuminated by a strong UV radiation field like in the galaxies in 
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our sample, HNCO and CH3OH are l argely photodissocia t ed an d only HCO, HOC + and 



C0 + should be observed like in M82 (IMartm et al. 



2006a 



20091 ). Then, not very massive 



molecular clouds and widely translucent to the UV radiation should dominate in M82. 
On the other hand, for galaxies with massive molecular clouds (large visual extinction) or 
low UV radiation fields, HNCO and CH3OH are well shielded and the abundance ratio 
of HNCO/CS will reach its maximum value. Considering that M83 and IC342 represent 
the stage of galaxies with an extremely low PDR component, the lower HNCO/CS ratios 
measured for NGC253, NGC4945 and NGC 1068 indicate that the PDR component must 
be substantial, as observed in other PDR tracers. 

For NGC 253 we find that the PDR component is similar to that in M82, but the total 
column density of dense gas is a factor of 2-3 larger in NGC 253 than in M82 from the low 
HNCO and CH 3 OH abundance in the latter. Considering the PDR column densities in both 
galaxies are similar, this component should represent about 1/3 — 1/2 of the total molecular 
column density in NGC253. This is roughly consistent with the decre ase by a factor of 



2 - 3 of the HNCO/CS ratio as compared with that of M83 or IC342 (IMartm et al.ll2009l ). 
This suggests that the molecular clouds properties in M82 and NGC 253 must be quite 
different in terms of the total molecular column density, which implies that the sizes or 
the densities are di fferent, or a c ombin atio n of both. Using t he atomic fine structure and 



CO emission lines, ICarral et al 



(119941 ) and 



Lord et all (119961 ) have also proposed that the 



clouds in M82 and NGC 253 are quite different. The clouds in NGC 253 are slightly smaller 
than those in M82, but with masses a factor of 15 larger than for M82. The NGC 253 
average cloud column densities are therefore factor of 20 larger than in M82. Though the 
total column densities of the M82 clouds inferred from the atomic fine structure lines are a 
factor of 5 larger than those predicted from the PDRs tracers and the HNCO abundance 
in this galaxy, similar constrains are derived both the molecular and the atomic tracers. 
Therefore, the clouds in NGC 253 are more massive than in M82. 
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We can even make a very rough estimate of the average properties and structure of 
the molecular clouds in NGC 253 and M 82 by combining the complementary information 



obtained from the 



from ( Martin et al. 



PR tracers presented in this paper and the HNCO column densities 



2009|). While HCO, HOC+ and CO + mainly trace the PDR region 
up to A v = 4 — 5 (i.e. H2 column densities of 5 x 10 21 cm -2 ), the HNCO emission only 
arises from the well shielded core (A v > 7) of the molecular clouds. The HCO, HOC + and 
CO + column densities in NGC 253 and M82 indicate similar averaged column densities 
in the PDR envelopes of the molecular clouds in both galaxies. The big difference in the 
molecular cloud structure in both galaxies is in the size (H2 column density) of the well 
shielded cores of the molecular clouds. In the case of M82, where HNCO has not been 
detected, we can set an upper limit to the HNCO column density of 7 x 10 12 cm~ 2 . This 
translates to a upper limit to the core H 2 column density of < 3 x 10 20 cm -2 for the HNCO 
fractional abundanc e of 2 x 10~ 8 derived for the well shielded clouds in the galactic center 



fjMartm et al. 



20081 ) . The averaged shielded cloud cores in M 82 are smaller by more than 
one order of magnitude than the PDR envelope. In the case of NGC 253, the H2 column 
densities of the shielded cloud cores is 10 22 cm -2 , a factor of 2 larger than PDR envelope. 
Assuming a similar averaged density distribution in the molecular clouds in both galaxies, 
the clouds in NGC 253 would be a factor 2 — 3 larger than in M82. 



4.4. The contribution X-ray induced chemistry 



Fuente et al. 



(120061 ) failed to repr oduce the large CO 



The PDR model presented by 
colum n density of a few 10 13 cm -3 observed towards M82. This lead 
(120071) to explore th e poss ibility of an enhanced X-ray induced chemistry in this galaxy. 



Spaans fc Meijerink 



Spaans Sz Meijerinkl (120071 ) concluded that such high formation of CO + can only be 
explained by X-ray irradiated molecular gas with densities of 10 3 — 10 5 cm~ 3 . Although the 
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X-ray luminosity of NGC 253 is a factor of 2 — 4 below th at of M82, both g alaxies have 



a significant X-ray emission in the range of ~ 10 40 ergs 1 (ICappi et al 



1999|). Similar to 



M82, the NGC 253 total C0+ column density is (3.6 ±1.1) x 10 13 cm" 2 . Moreover, both 
show a similar HCO+ /CO + ratio of - 30 - 40. 

The models presented in this paper are able to produce such column densities for visual 
extinctions of A v ~ 3 — 5 for model A, and A v ~ 0.5 — 1 for model B. These models have been 
calculated with a radiation field, a cosmic ray flux and a H2 density sm aller by a factor of 2 , 



Fuente et al. 



(|2006|). 



40 and 4, respectively, with respect to those assumed in the models of 
Furthermore, we are able to reproduce the abundances and abundance rat i os me asured for 
all the other observed species presented in this paper. 



van der Tak et al. 



(120081 ) showed 



how the measured abundance of H3O 4 " in M82 can be both produced by PDR with a high 
cosmic-ray ionization or by an XDR. Indeed, an increase of cosmic ray ionization rate in 
PDR models may be qualitatively used to simulate XDR-like environments. Our models, 
however, do not use particularly high cosmic ray fluxes (a factor of 5 higher than standard). 
Thus, though the X-ray irradiation is substantial in SB galaxies, the PDR models presented 
in this paper can reproduce the molecular abundances observed towards the brightest 
prototypes, M82 and NGC 253. Moreover, no significant changes in the abundances of 
HOC + and HCO are found towards the nuclear AGN of NGC 1068 where X-ray radiation is 
significantly more important than in SB nuclei, as shown in Sect. 14.1.21 Unfortunately, no 
CO + observation has been reported towards this Seyfert 2 nucleus. 



5. Conclusions: The pervading UV field in evolved starburst 

The comparison of model predictions with the observations presented show that 
the abundance of the species observed in this work towards NGC 253, namely HCO + , 
CO + , and HCO, are most efficiently formed in the outer region of the molecular clouds 
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where the gas is highly irradiated by the incident UV photons from massive stars. The 
high molecular abundances derived for these species in NGC253 suggest that the PDR 
component in this galaxy is similar to that found in M82, claimed to be the prototype 
of extragalactic PDR. The abundance ratios found for this limited sample of galaxies are 
of the same order as those observed towards galactic PDRs, which stress the importance 
of photo-dominated chemistry in galaxy nuclei. Large amounts of molecular material are 
affected by photodissociation not only in NGC253, but also towards the star forming 
regions around the Seyfert 2 nuclei in NGC4945 and NGC1068. This is consistent with 
the HNCO / CS ratio in these galaxies which suggest that a fraction of HNCO has been 
photodissociated in PDRs. The combination of the observations of HCO, HOC + and CO + 
with that of HNCO seems to confirm that their abundances reflect the evolutionary stage of 
the starbursts in these galaxies. Although photodissociation is the most likely scenario for 
the enhancement of the observed reactive ion in starburst environments, X-ray dominated 



AGNs in circunnuclear disk of NGC 1068 fHOC+ 
luminous infrared galaxy Arp 220 (H 3 



(HOC+ 


Usero et al. 


2004 


) and towards the ultra 


van der Tak et al. 


2008 


)• 



Therefore, M82 is still outstanding not only as a PDR dominated galaxy, 
but by the underabundance o f complex mo l ecules such as CH 3 OH, HNCO or SiO 



( IMauersberger fc Henkel 



1993 



Martin et al 



2006a 



bj), evidence for the lack of large 



amounts of dense molecular material which would pot entially fuel its nuc lear starburst 



combination with the HNCO abundances (IMartm et al. 



as compared to other starburst galaxies li ke NGC 253 ([Martin et al. 



20091 ). Our data in 



2009) indicate that the molecular 



clouds in M82 are different from those in NGC253. Although having a similar overall 
PDR component, the clouds in NGC 253 have to be more massive and have larger column 
densities those in M 82. 
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Fig. I.— IRAM 30m observations of the emission of HCO+, C 18 and HOC+ J = 1 - 
transitions with the fitted line profiles superimposed. C 18 emission has been multiplied by 
a factor of 4 for comparison with the HCO + profile. The spectral resolution are the original 
1 MHz (~ 3 kms -1 ) and smoothed to 4MHz (~ 13 kms -1 ), respectively. 
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Fig. 2. — HCO profile in the same window as SiO and H 13 CO + observed with the IRAM 
30 m.. On top it is shown the triple Gaussian profile fitted to the transitions. Dashed lines 
indicate the position of each hyperfine structure transition of HCO, where only the brighter 
one has been fitted. The theoretical synthetic spectrum of HCO is shown with a dotted line 
(see Sect. [2] for details). Velocity resolution has been degraded to ~31kms _1 . 
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Fig. 3. — JCMT observations of the CO + emission blended with transitions of CH3OH J = 
5 — 4 and observed in the same window as HC3N J = 25 — 24. This is the first detection 
of the 13 C isotopologue of methanol. The overall spectral fitting to all the lines is shown 
in grey. The contribution of the CO + emission is shown with dashed line. The position of 
the three brighter transitions of the J = 5 — 4 group of 13 CH 3 OH are shown indicated with 
vertical dashed lines. See Sect. [2] for details on the fitting to the spectra. Velocity resolution 
has been degraded to ~30kms _1 . 
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Fig. 4. — Theoretical predictions for the fractional abundances relative to H 2 ( Upper panels) 
and abundance ratios (Lower Panels) for the observed species as derived from the two dif- 
ferent PDR models: a pure gas-phase model A (Left panels) and a coupled dense core-PDR 
model B (Right panels). Details are given in Section I4~2l The vertical position of the key for 
each molecule and ratio shown only in the plots for model A correspond to the actual derived 
parameters from the observations. Additionally, the fractional abundances of CH3OH and 
HNCO are shown for model B. 
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Table 1: Parameters derived from the observed line profiles. 



Transition 




/ T MB dv 


^LSR 


A«i/2 


Tub 






(Kkms" 1 ) 


(kms -1 ) 


(kms^ 1 ) 


(mK) 


C 18 f - 




3.1 ±0.6 


183 ±13 


100 ±9 


116.75 






6.0 ±0.8 


283 ±7 


100 ± 9 a 


223.48 


HCO+ 1 - 




21.98 ± 1.1 


177.9 ±0.4 


118.1 ±0.3 


174.9 






35.79 ±1.4 


289.0 ±0.3 


118.1 ±0.3 a 


284.8 


HOC+ 1 - 




0.8 ±0.2 


170 ±10 


100 ±20 


7.6 






1.0 ±0.2 


282 ±9 


100 ±20 


9.2 


HOC+3 - 2 




< 0.8 b 






< 4.1 


HCO l ,i - ,o 




0.41 ±0.08 


183 ± 14 


102 ± 4 a 


3.8 






0.43 ±0.08 


297 ± 13 


102 ± 4 a 


3.9 


H 13 CO+ f - 




1.26 ±0.09 


176 ±6 


102 ± 4 a 


11.6 






1.36 ±0.10 


285 ±5 


102 ± 4 a 


12.5 


SiO 2 - 1 




1.52 ±0.11 


182 ±5 


102 ± 4 a 


13.9 






1.59 ±0.11 


292 ±5 


102 ± 4 a 


14.6 


HC 3 N J = 5 - 4 




1.05 ±0.15 


191 ±7 


95 ± 18 


10.5 


CO+ 5/2 - 3/2F = 


2 - f 


0.30 ±0.10 


191 c 


95 c 


3.0 


CO+ 3/2- 1/2F = 


2 - f 


0.17 ±0.10 


191 c 


95 c 


1.6 


13 CH 3 OH 5 ,5 - 4 .4 




0.10 ±0.07 


191 c 


95 c 


1.0 


13 CH 3 OH 5_i,5 -4_ 


-1.4 


0.16 ±0.07 


191 c 


95 c 


1.6 



11 Linewidths forced to have the same value in the Gaussian fit. 

3tr upper limit assuming a 200 kms -1 lincwidth. 
c Parameters forced to equal those derived from HC3N Gaussian fit. 
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Table 2: Derived fractional abundances and H 13 CO + ratios for each velocity component 



Molecule 


N 




H 13 CO+/X 




(xlO 13 cm" 2 ) 


(10~ 10 ) 




H 13 CO+ 


1.6 ±0.5 


4.8 ± 1.0 


1 




1.7 ±0.6 


2.7 ±0.4 


1 


HOC+ 


0.8 ±0.3 


2.4 ±0.7 


2.0 ±0.8 




1.1 ±0.4 


1.7 ±0.4 


1.6 ±0.4 


HCO 


12.3 ±5.8 


37 ± 10 


0.13 ±0.04 




12.9 ±6.1 


20 ±4 


0.14 ±0.03 


CO+ 


1.7 ±0.8 


5±2 


0.9 ±0.4 



a With N(H 2 ) = 3.3 ± 1.2 x 10 22 cm 2 and 6.3 ±2.1 x 10 22 cm 2 for each velocity component, respectively, as 
derived from C 18 with 16 0/ 18 = 150 (jrlarrison et al.lll999l ) and a CO/H 2 = 1CT 4 . 
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: Abundance ratios of HCO + vs 


HOC+ and HCO. 




Source 


l~T T / \ _L_ 1 / |"T T / \ /~~~t -4- 1 

[HCO + ]/[HOC + ] 


[HCO+]/[HCO] 


[HCO+]/[CO+] 


JNGC 253 


80 ±30 


5.2 ± 1.8 


38 ± 15 




63 ± 17 


5.4 ±1.3 




M82 


60±28 a 


9.6±2.8 b 


32 ± 16 a 


NGC 1068 


128 ±28 c 


3.2 ± 1.2 d 




NGC 4945 




2.4 ± 1.2 C 






GC prototypical PDRs 




Horsehead 


75 - 200 f 


1.1 g 


> 1800 f 


Orion Bar 


< 166 - 270 h ' 1 


2.4 J 


< 83 - 140 


NGC 7023 


50 - 120 1 


31 J 


14 1 


S140 


4110 k 


3.5- > 62 U 


7420 k 


NGC 2023 


940 k 


1.1 j 


740 k 


M17SW 


2262 h 




250 m 



a Derived from single dish data (jMauersberger fc Henkellll991t iFuente et al 
b Average ratio from the interferometric maps by ( Garcia-Burillo et al 



: Average over the whole line profile in the CND position (jUsero et al 



200^ 



2003 1. See Sect. O for details. 



20ol . 



d Aver age value over the three positions in the circunnuclear starburst ring with HCO detections (jUsero et al 



20041) 



Wang et al 



e From 

f Goicoechea et al 



( 20041) 



(|2009h 



Gerin et al 



( 20091) 



Apponi et al 



( 19991) 



Orion Bar ionization front and PDR-peak in NGC 7023 



Fuente et al 



( 20031) 



Schilke et al. 


(200 


1) 


Savage & Ziurvs ( 


2004) 


Schenewerk et al. 


(1988) 


Storzer et al. 


(199 


5) 



